function [act logact] = gaussBayes(m, b, W, v, X)


[N D] = size(X);

centered_X = X - ones(N, 1)*m;
exp_dist2 = D/b + v*sum(centered_X*W.*centered_X, 2);
exp_logdet = wishart_exp_logdet(W, v);

logact =  -0.5*D*log(2*pi) + 0.5*exp_logdet - 0.5*exp_dist2;
act = exp(logact);
